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Introduction 

Transfer  of  heat  into  the  soil  beneath  a  wildland 
fire  causes  a  large  number  of  onsite  fire  effects 
(Hungerford  1990)  that  include  mortality  of  or  in- 
jury to  subsurface  parts  of  surface  plants  or  whole 
organisms  that  live  within  the  soil,  thermolysis  of 
organic  matter,  oxidation  or  volatilization  of  free 
mineral  components,  and  other  physicochemical 
changes  that  can  be  as  extreme  as  fusion  of  mineral 
soil  constituents.  To  predict  the  nature  and  extent  of 
such  effects  it  is  necessary  to  know  the  temperature 
history  as  a  function  of  depth  within  the  soil  beneath 
the  fire. 

To  predict  the  temperature  history  profile  in  a  soil 
under  a  fire,  we  need  to  know  the  history  of  heat  in- 
put to  the  soil  surface  or  the  surface  temperature 
history  and  to  model  the  response  of  the  soil  to  this 
stimulus.  In  this  paper  we  focus  on  the  challenge  of 
modeling  the  response  of  the  soil  to  boundary  heat- 
ing or  temperature  history. 

Modeling  the  soil's  thermal  response  would  be 
relatively  simple  if  not  for  moisture  and  its  move- 
ment because  of  spatially  nonuniform  heating  (Morse 
and  Feshbach  1953;  Steward  and  others  1990).  The 
very  large  heat  of  vaporization  of  water  and  the  ease 
with  which  water  vapor  can  move  through  a  porous 
soil  matrix  cause  the  transport  of  heat  to  be  domi- 
nated by  water  vapor  movement  for  temperatures 
near  the  boiling  point  of  water  (Aston  and  Gill  1976; 
Cheng  1978;  Eckert  and  Faghri  1980;  Luikov  1975). 
The  relative  immobility  of  water  as  a  liquid  in  soil, 
and  the  high  value  of  hquid  water's  thermal  conduc- 
tivity relative  to  that  of  air  in  the  soil  pores,  make 
the  thermal  conductivity  of  the  bulk  soil  medium 
sensitive  to  its  liquid  water  content  when  liquid  wa- 
ter occupies  a  significant  volume  fraction  within  the 
medium  (Campbell  1985;  deVries  1958;  Jury  1973; 
Peter  1992;  Philip  and  deVries  1957). 

All  this  complicates  the  accurate  modeling  of  heat 
transport  within  soils  by  making  it  necessary  simul- 
taneously to  model  the  movement  and  concentration 


of  water  in  liquid  and  vapor  phases  as  well  as  the 
transport  of  heat  by  thermal  conduction.  Further- 
more, a  complete  model  must  include  variations  with 
moisture  content  and  temperature  of  the  param- 
eters that  control  the  transport  of  heat,  liquid  wa- 
ter, and  water  vapor.  The  dependence  of  rate-con- 
trolling parameters  (such  as  thermal  conductivity, 
vapor  dififusivity,  and  hydraulic  conductivity)  upon 
the  dependent  variables  means  that  the  equations 
are  mathematically  nonlinear.  This  nonlinearity, 
makes  it  not  valid  to  solve  a  series  of  simple  prob- 
lems and  construct  more  complex  cases  by  linear  su- 
perposition of  the  "simple"  solutions.  Nonlinear  sys- 
tems of  equations  do  not  generally  admit  to  analytical 
solutions.  Thus,  we  need  to  resort  to  numerical  inte- 
gration of  the  model's  partial  differential  equations, 
introducing  a  whole  new  realm  of  potential  pitfalls. 
Despite  these  complexities,  a  fairly  rich  literature 
on  the  topic  exists. 

Heat  transport  and  moisture  movement  within 
soils  have  been  studied  in  soil  science  and  micro- 
climatology  (Al  Nakshabandi  and  Kohnke  1965; 
Bhuiyan  and  others  1971;  Campbell  1985;  deVries 
1958;  deVries  and  Afgan  1975;  Philip  and  deVries 
1957).  There  is  a  wealth  of  theoretical  and  empirical 
knowledge  about  this  phenomenology  and  the  mod- 
eling of  it.  But  the  models  and  relationships  are  fi-e- 
quently  restricted  to  the  temperature  regimes,  heat- 
ing rates,  and  spatial  gradients  of  temperature  and 
moisture  that  occur  naturally  in  diurnal  cycles  of 
heating  and  cooling  (and  wetting  and  drying)  and 
as  a  result  of  precipitation,  drought,  and  irrigation. 

Only  recently  (Hungerford  1990)  have  investiga- 
tions been  extended  to  the  temperatures,  heating 
rates,  dryness  states,  and  gradient  magnitudes  that 
characterize  the  fire  environment  (Aston  and  GiU 
1976;  Campbell  and  others  1992,  1993, 1994,  1995; 
Jury  1973;  Peter  1992;  Schroeder  1974;  Steward 
and  others  1990).  However,  none  of  these  authors 
have  drawn  heavily  upon  parallel  literature  fi-om 
the  disciplines  of  engineering  and  geophysics. 
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In  those  research  communities,  analytical,  numerical, 
and  experimental  investigations  of  coupled  heat  and 
mass  transport  in  porous  media  have  been  carried 
out  for  more  than  a  century  (Catton  1992;  Kaviany 
1991;  Rust  and  Roberts  1990;  Tang  and  Bau  1992). 
Engineers  view  this  as  a  special  case  in  the  general 
field  of  heat  transfer  (Bejan  1984).  Half  of  this  re- 
cent volume  (Kaviany  1991)  is  devoted  to  analysis 
of  heat  and  mass  transport  in  the  instance  that  the 
interstitial  fluid  exists  in  two  phases,  which  is  the 
case  of  interest  in  fire-heated  soils. 

This  paper  documents  a  research  joint  venture. 
The  objectives: 

1.  Survey  the  literature  across  several  disciplines 
that  address  heat  and  mass  transport  in  porous  me- 
dia, assessing  any  differences  in  approach  to  model- 
ing these  processes  among  the  different  disciplines. 

2.  Document  test  exercises  of  two  models  that  are 
potentially  useful  for  assessing  effects  of  wildfires. 

3.  Identify  opportunities  for  advancing  the  state  of 
the  art  in  modeling  fire-driven  transport  of  heat  and 
moisture  in  soils. 


Survey  of  Literature 

A  literature  survey  was  undertaken  through  the 
technical  libraries  at  Montana  State  University  and 
at  Intermountain  Fire  Sciences  Laboratory,  as  well 
as  the  authors'  personal  literature  collections  and 
preprints  fi"om  technical  meetings.  We  found  de- 
scriptions of  computer  simulation  models  for  heat 
and  mass  transport  in  porous  media  in  the  litera- 
ture of  soil  science  and  its  derivative  subdisciplines, 
and  in  publications  by  engineers  and  by  geophysicists. 

Models  with  Soil  Science  Ancestry 

Hungerford  (1990)  offers  a  definitive  statement  of 
the  motivation  for  modeling  heat  and  moisture  trans- 
port in  soils  for  the  purpose  of  predicting  the  effects 
of  wildland  fires  and  provides  a  comprehensive  sur- 
vey of  empirical  and  interpretive  field  studies  of  the 
effects  of  fires  in  forest  and  rangeland  settings.  He 
provides  an  exhaustive  list  of  the  mathematical  mod- 
els that  had  been  put  forward  by  mid- 1990  for  the 
prediction  of  heat  and  moisture  transport  in  soils 
heated  by  fires. 

The  pioneering  works  by  J.  R.  Philip  and  D.  A. 
deVries  (deVries  1958;  deVries  1975;  Phihp  1975; 
Philip  and  deVries  1957)  explored  and  analyzed  the 
movement  of  moisture  and  the  transfer  of  heat  in 
soils  by  starting  at  the  size  scale  of  the  soil  particles. 
They  sought  to  connect  soil  descriptive  parameters 
such  as  the  relative  abundances  of  mineral  types  in 
the  soil,  particle  shapes  and  size  distributions,  bulk 


density,  soil  environmental  parameters  such  as  tem- 
perature and  liquid  moisture  volume  fraction,  and 
the  constitutive  parameters  describing  macroscale 
phenomenology;  such  as  hydraulic  conductivity,  ef- 
fective vapor  dififusivity,  and  apparent  thermal  con- 
ductivity. This  approach  has  had  considerable  suc- 
cess and  is  the  current  paradigm  of  the  soil  science 
commvmity  for  modeling  heat  and  moisture  move- 
ment in  response  to  diurnal  fluctuations  of  the 
forces  driving  them. 

deVries  Formulation— In  a  siu-vey  paper,  deVries 
(1958)  summarized  and  generalized  simultaneous 
differential  equations  for  the  transfer  of  heat  and 
of  moisture  presented  the  year  before  (Philip  and 
deVries  1957)  but  is  available  in  restricted-circula- 
tion publications  and  in  a  fragmented  form. 

We  found  no  derivations  of  the  moisture  and  heat 
transport  equations  in  either  of  these  citations.  The 
conservation  of  mass  laws  for  water  vapor  and  for 
liquid  water  are  used  with  models  for  fluxes  embed- 
ded when  they  are  stated,  and  an  energy  conserva- 
tion equation  is  assembled  around  the  Fourier  tran- 
sient heat  conduction  equation.  The  later  paper  is 
summarized  briefly  here,  emphasizing  a  few  confus- 
ing or  inconsistent  aspects  of  this  important  work. 

Phenomenological  equations  relate  moisture 
fluxes  in  liquid  and  vapor  phases  to  gradients  in 
temperature  and  moistxire  content.  These  relation- 
ships would  supposedly  arise  from  the  separation  of 
the  equation  expressing  conservation  of  momentum 
for  the  medium  considered  as  a  continuum  and  then 
decomposed  into  its  parts:  soil  particles,  liquid  wa- 
ter, water  vapor,  and  interstitial  air.  This  would  be 
the  paradigm  in  an  engineering  approach. 

Two  critical  assumptions  appear  early  in  the  paper: 

1.  The  moist  soil  medium  is  microscopically  iso- 
thermal (all  components  in  the  near  neighborhood 
of  a  location  are  assigned  the  same  temperature)  yet 
large  temperature  gradients  can  exist  across  pore 
spaces  between  solid  particles.  deVries  does  not  dis- 
cuss the  contradiction  in  this  two-part  assumption. 

2.  Water  vapor  and  liquid  water  are  in  local  vapor 
pressure  equilibrium.  This  assumption  allows  the 
combining  of  liquid  water  content  and  water  vapor 
content  to  express  a  water  content  distribution  un- 
ambiguously, in  absence  of  hysteresis,  as  the  appor- 
tioning between  vapor  and  liquid,  then  depends 
uniquely  upon  porosity  and  local  temperature. 

While  neither  vapor  buoyancy  nor  liquid  viscosity 
are  found  explicitly  in  the  parameters  given,  the  bi- 
nary molecular  diffusion  coefficient  of  water  vapor 
in  air  appears  in  proportionality  constants  for  both 
thermally  induced  movement  of  vapor  and  movement 
of  vapor  in  response  to  gradients  in  water  potential. 
The  rate  of  evaporation,  expressed  as  rate  of  change 
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of  mass  of  water  (measured  as  equivalent  volume  of 
liquid  water)  per  unit  volume  per  unit  time,  appears 
as  a  source  term  in  the  primitive  foirn  of  the  vapor 
motion  equation  and  as  a  sink  term  in  the  liquid 
motion  equation.  Instead  of  relating  the  rate  of  evap- 
oration to  the  net  absorption  rate  of  energy,  after 
accounting  for  sensible  heat  accumulation,  the  term 
is  eliminated  by  summing  the  equation  for  water  va- 
por conservation  and  the  equation  for  liquid  water 
conservation  to  form  an  equation  for  global  conser- 
vation of  water.  This  step  illuminates  a  difference 
between  the  soil  science  approach  and  that  of  the 
investigator  of  the  effects  of  fire. 

The  flux  of  heat  is  related  to  the  transport  of  sen- 
sible and  latent  heat  by  water  vapor,  transport  of 
sensible  heat  by  liquid  water  movement,  and  heat 
conduction  through  the  continuum.  At  first,  conduc- 
tion through  the  continuum  is  modeled  by  a  "hypo- 
thetical" thermal  conductivity,  X*,  understood  to 
characterize  soil  in  which  the  moisture  distribution 
is  fixed  (deVries  1958,  eq.  10).  Later,  this  symbol  is 
replaced  by  the  collection: 

X,  -  L  pi 

where  X  is  identified  as  the  "apparent  conductivity 
of  the  medium  including  the  effect  of  vapor  distilla- 
tion due  to  temperature  gradient,"  L  is  the  latent 
heat  of  vaporization  of  water,  pi  the  density  of  liquid 
water,  and  Dtv  the  thermal  diffusivity  of  water  vapor 
within  the  soil  medium.  The  reader  is  cautioned 
that  the  two  quantities  are  not  equal  because  the 
heat  fluxes  in  different  parts  of  the  medium  are  not 
additive. 

An  investigation  of  steady  state  heat  conduction 
concludes  the  paper  and  illustrates  the  model's  non- 
linear character.  A  plot  of  isothermal  moisture  con- 
ductivity (velocity  of  liquid  water  movement  per  imit 
gradient  of  water  potential,  with  potential  expressed 
as  hydraulic  head)  and  isothermal  moisture  diffu- 
sivities  (as  vapor  and  liquid)  in  relation  to  gradients 
in  moisture  content  show  strong  sensitivity  to  water 
content.  Another  plot  similarly  shows  isothermal 
liquid  water  potential  as  a  function  of  liquid  mois- 
ture content  and  moisture  diffusivities  in  response 
to  thermal  gradients  as  strongly  dependent  upon 
liquid  moisture  content. 

Using  the  mathematical  relationships  illustrated 
in  these  two  figures,  deVries,  found  steady  state 
equations  for  temperature  gradient  and  liquid  mois- 
ture gradient  by  setting  to  zero  the  time  dependent 
terms  (local  temperature  rise  rate  and  moisture  ac- 
cumulation rate)  and  the  moisture  flux.  He  then 
solves  these  steady  state  equations  for  the  steady 
state  spatial  temperature  and  moisture  gradients 
as  a  function  of  liquid  moisture  content.  The  final 
two  figures  in  the  paper  display  the  results  of  these 
manipulations  for  a  light  clay  and  for  a  medium 


sand  son,  at  20  °C  with  a  constant  horizontal  heat 
flux  of  10"^  cal/s  cm^  (about  42  W/m^).  The  moisture 
and  temperature  distributions  in  space  are  not 
given  but  can  be  inferred  qualitatively  fi-om  these 
graphs  to  be  approximately  as  follows:  For  the  light 
clay  soil,  volumetric  moisture  content  rises  exponen- 
tially with  distance  fi*om  the  plane  of  heat  applica- 
tion, with  a  characteristic  length  of  about  10  cm. 
When  the  liquid  moisture  content  reaches  about  0.05, 
the  slope  of  this  curve  increases  by  a  factor  of  about 
15;  the  moisture  content  rises  abruptly  and  roiinds 
off  to  its  asymptotic  value  of  about  0.12.  The  physi- 
cal processes  that  enforce  such  a  variation  of  mois- 
ture content  with  distance  are  not  extensively  dis- 
cussed. For  a  medium  sand,  the  moisture  content 
varies  initially  with  distance  fi-om  the  plane  of  heat 
apphcation  again  as  an  exponential,  but  with  a  char- 
acteristic length  nearer  1  mm.  Upon  reaching  a  hq- 
uid  moisture  content  of  about  0.005,  the  cvirve  flat- 
tens abruptly,  the  distance  scale  of  the  exponential 
growing  to  about  80  cm.  This  profile  would  resemble 
a  slightly  round-shouldered  step  function. 

The  parameter  values  used  by  deVries  are  repre- 
sentative of  the  concerns  of  soil  scientists  in  quanti- 
fying phenomenology  associated  with  irrigating,  the 
drying  of  soils  heated  by  sunshine,  and  evaporation 
of  soil  moisture  into  a  dry  atmosphere.  The  heating 
rates  and  temperature  levels  of  interest  in  soils  un- 
der fires  are,  of  course,  substantially  different.  Some 
of  the  features  of  the  model  presented  by  deVries 
are  confusing  or  surprising  probably  due  to  the  dif- 
ferences in  technical  backgrounds  and  in  motivation 
for  addressing  this  complex  phenomenology.  Had  an 
engineer  specializing  in  fluid  mechanics  or  in  heat 
transfer  described  the  processes,  the  terms  used  and 
the  emphasis  would  be  much  different.  An  engineer 
would  have  selected  different  ranges  of  parameters, 
if  the  perturbation  of  soil  temperature  and  moisture 
profiles  were  to  be  inferred  due  to  heating  by  a  fire 
in  the  surface  mantle. 

Most  modem  engineering  audiences  would  find 
the  deVries  (1958)  formulation  deficient  in  at  least 
two  respects:  First,  the  model  contains  no  equation 
for  the  conservation  of  momentum.  Because  there  is 
no  explicit  expression  of  a  force  balance,  one  is  left 
to  ponder:  Does  the  water  vapor  go  up  or  down  when 
the  soil  is  heated  on  the  upper  surface?  The  vapor 
tends  to  move  from  high  temperature  toward  low 
temperature,  but  also  to  move  from  high  toward  low 
moisture  concentration.  The  mechanical  forces  are 
not  identified  except  in  limiting  forms,  nor  does  the 
acceleration  of  gravity  appear  in  the  model.  Thus, 
the  phenomenon  that  nonuniform  pressure  acceler- 
ates a  fluid  toward  lower  pressure  is  omitted,  so  wa- 
ter vapor  has  no  buoyancy  or  tendency  to  rise.  Con- 
ditions under  which  the  neglect  of  buoyancy  is  valid 
need  to  be  established. 
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Second,  the  energy  conservation  equation  is  not 
derived  from  the  energy  transport  equation  of  con- 
tinuum mechanics,  so  the  transfer  of  energy  by 
mass  movement  is  "added  on"  to  the  Fourier  heat 
conduction  equation.  Such  a  derivation  taxes  com- 
prehension because  the  heat  conduction  equation 
employs  an  apparent  thermal  conductivity  that  in- 
cludes the  transport  of  heat  by  vapor  movement. 
This  contribution  is  then  subtracted  from  conductiv- 
ity measurement  data  so  the  effect  can  be  added  in 
to  the  vapor  movement  model.  This  formulation  does 
not  readily  permit  differences  in  temperature  between 
liquid,  vapor,  and  mineral  constituents,  which  would 
be  relatively  simple  to  accommodate  in  a  model  de- 
rived from  the  perspective  of  continuum  mechanics. 

Although  these  aspects  of  deVries'  (1958)  formula- 
tion are  not  discussed  in  later  literature,  this  water- 
shed work  gave  us  a  model  for  the  coupled  transport 
processes  that  has  persisted  largely  unchanged,  ex- 
planatory power  must  be  acknowledged  despite  any 
uncertainties  raised  here. 

Fire  Heating  Model  of  Aston  and  Gill— Another 
early  and  important  model  is  that  of  Aston  and  Gill 
(1976),  who  addressed  the  problem  for  conditions  of 
heating  by  a  surface  fire.  These  authors  address  the 
modeling  of  heat  and  moisture  transport  in  soils  in 
the  vertical  direction  only,  but  generalization  of 
their  equations  should  be  straightforward. 

Compared  to  deVries'  complicated  formulation, 
that  of  Aston  and  Gill  embodies  simplicity  even 
though  it  rests  upon  three  coupled  diffusion  equa- 
tions, each  with  a  source  term.  Equations  are  pos- 
ited as  phenomenological  descriptions  of  sensible 
heat  transport,  liquid  water  transport,  and  water 
vapor  transport.  Not  emphasized  is  the  assumption 
that  all  components  share  a  local  temperature. 

The  transport  of  sensible  heat  is  modeled  without 
mass  transport  except  that  the  rate  of  condensation 
of  water  vapor,  multiplied  by  the  latent  heat  of  va- 
porization, appears  as  a  source  term  in  the  tran- 
sient Fourier  equation  for  heat  conduction.  In  this 
equation,  the  specific  heat  capacity  per  unit  volume 
is  a  function  of  the  porosity  and  liquid  water  content 
of  the  soil,  and  the  conductivity  of  the  medium  as  a 
continuum  is  related  empirically  to  the  same  pa- 
rameters, following  Al  Nakshabandi  and  Kohnke 
(1965).  Transport  of  sensible  heat  by  moisture  (as 
liquid  or  vapor  or  both)  is  incorporated  through  use 
of  an  empirical  thermal  conductivity  (deVries'  "ap- 
parent" thermal  conductivity),  for  which  measure- 
ments usually  exist  only  at  low  temperatures.  The 
authors  used  a  linear  function  of  temperature  as  an 
empirical  correction  factor  of  the  thermal  conductiv- 
ity, allowing  them  to  match  measured  histories  of 
temperature  at  various  depths  beneath  a  surface 
exposed  to  radiant  heating. 


The  movement  of  liquid  water  is  modeled  by  the 
continuity  equation,  the  speed  of  liquid  water  mo- 
tion being  given  by  the  product  of  the  gradient  of 
hydraulic  potential  and  the  hydraulic  conductivity. 
A  source  term  is  included,  being  the  local  rate  of 
condensation  of  water  vapor.  We  believe  that  the 
"hydraiilic  potential"  described  here  is  rather  the 
water  (matric)  potential  expressed  as  equivalent 
water  depth  at  standard  gravity. 

The  movement  of  water  vapor  is  modeled  by  a  con- 
tinuity equation,  with  local  velocity  of  water  vapor 
modeled  as  being  equal  to  the  product  of  the  gradi- 
ent of  vapor  pressure  and  "vapor  conductivity"  in 
analogy  to  hydraulic  conductivity  for  the  liquid 
phase.  This  parameter  is  proportional  to  the  binary 
molecular  diffusion  coefficient  and  includes  an  em- 
pirical correction  factor  to  be  chosen  by  the  model 
user. 

The  vapor  pressure,  e,  is  derived  from  the  tem- 
perature-dependent saturation  vapor  pressure,  e^, 
scaled  by  a  temperature  sensitive  exponential  factor: 

ele^  =  e^^i-^MJRT) 

where 

<D    =  "hydraulic  potential"  (must  be  converted 
to  erg/gm) 
-  molecular  weight  of  water,  gm/mol 
R    =  universal  gas  constant,  erg/mol  K 
T    =  absolute  temperature,  K 

Note  the  conflict  in  dimensions  of  the  parameters 
stemming  from  the  fact  that  the  unit  of  the  "hydrau- 
lic potential"  is  given  as  length.  If  the  acceleration  of 
gravity  is  added  as  a  factor  in  the  numerator  of  the 
exponent,  it  becomes  dimensionless.  This  is  consis- 
tent with  the  formulation  of  Bhuiyan  and  others 
(1971),  cited  by  Aston  and  Gill  as  the  source  for  the 
liquid  water  motion  model  used.  Bhuiyan  and  others 
use  the  same  symbol,  O,  but  identify  it  as  the  sum  of 
the  "pressure  potential"  and  the  "gravitational  po- 
tential," the  latter  having  the  units  of  acceleration 
times  length. 

The  authors  did  not  use  the  vapor  pressiire  formal- 
ism in  the  numerical  examples  they  presented.  In- 
stead, e/Cg  was  assumed  to  be  independent  of  tem- 
perature and  was  taken  from  an  empirical  equation 
relating  it  to  volumetric  water  content  for  a  particu- 
lar soil.  Again  a  source  term  was  included  in  the  liq- 
uid water  continuity  equation,  being  the  negative  of 
the  water  vapor  source  term. 

Boundary  conditions  used  for  the  computer  imple- 
mentation of  this  model  are  in  the  form  of  a  specified 
surface  temperature  history.  The  moisture  bound- 
ary condition  is  unspecified  at  the  heated  surface, 
but  temperature  and  moisture  content  are  held 
constant  at  a  "deep"  boundary,  where  their  gradi- 
ents vanish.  Numerical  results  generated  from  the 
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computer  code  of  Aston  and  Gill  (1976)  showed 
agreement  with  soil  temperatures  measured  vmder 
a  spreading  grass  fire  (Scotter  1970).  The  tempera- 
ture predictions  of  the  model  range  from  room  tem- 
perature to  over  400  °C,  moisture  contents  treated 
range  from  bone  dry  to  near  saturation,  and  simu- 
lated time  extended  well  beyond  200  minutes.  These 
results  are  impressive  considering  that  the  model 
equations  are  mathematically  "stiff"  (Press  and  oth- 
ers 1986),  and  especially  in  light  of  the  questionable 
features  noted  in  the  model  for  heat  transport.  Al- 
though this  model's  predictions  agreed  well  with  the 
data  of  Scotter  (1970),  the  model  has  not  performed 
well  when  applied  under  other  conditions,  and  it  ap- 
pears to  lack  generality. 

A  Model  Ignoring  Moisture  Movement— In 

using  a  model  of  heat  transport  in  soil,  one  common 
prediction  is  whether  certain  plant  components  are 
killed  by  the  thermal  impact  of  a  fire.  Most  plant  tis- 
sue (excepting  seed)  is  killed  once  its  temperature  is 
raised  to  60  °C  for  a  short  time  (Levitt  1980).  How 
far  below  the  surface  of  the  mineral  soil  at  a  fire  site 
is  this  maximum  temperature  achieved?  The  sim- 
plest model  for  predicting  that  maximal  depth  dur- 
ing a  spreading  fire  is  one  that  ignores  moisture 
movement  within  the  soil  and  treats  the  moist  me- 
dium as  an  inert  solid  with  constant  constitutive  pa- 
rameters. Modeling  the  fire  as  a  moving  line  heat 
source  at  the  soil  surface,  Richon  (1987)  found  that 
for  the  soil  surrogate  material  properties  explored, 
the  problem  could  be  simplified  by  idealizing  the 
analysis  to  a  one-dimensional  transient  situation. 
Steward  and  others  (1990)  extended  this  idealiza- 
tion to  include  cooling  of  the  heated  surface  using 
a  constEint  Newtonian  film  heat  transfer  coefficient. 
They  found  analytical  solutions  in  dimensionless 
form  for  four  different  profiles  of  heat  transfer  rate 
as  a  function  of  time,  with  Newtonian  cooling  occur- 
ring either  at  £ill  times  or  only  after  heating  had 
ceased. 

Soil  can  be  treated  as  a  medium  with  fixed  prop- 
erties and  the  fire  treated  as  a  transient  bovmdary 
heating  rate  (to  the  extent  that  such  approximations 
do  not  impose  intolerable  error),  a  treatment  such  as 
this  is  highly  desirable.  Using  dimensionless  vari- 
ables, aU  situations  of  an  entire  class  can  be  described 
by  a  single  equation. 

We  cannot  solely  rely  on  models  using  moist  soil 
as  a  medium  with  constant  thermophysical  proper- 
ties. However,  one  must  assume  that,  even  though 
all  the  assumptions  of  the  model  are  not  valid,  it 
may  have  adequate  predictive  power.  Nevertheless 
extensive  empirical  justification  is  not  yet  available 
for  this  model. 


CampbelFs  Model— Prominent  soil  scientist, 
Professor  Gaylon  Campbell  of  Washington  State 
University,  led  a  research  team  in  developing  a 
model  for  heat  and  moisture  transport  in  soils  appli- 
cable at  the  heating  rates  and  temperature  levels 
of  soils  under  wildland  fires.  This  model  (Campbell 
and  others  1992,  1995)  embodies  several  significant 
extensions  of  deVries'  formulation.  It  includes  equa- 
tions for  predicting  apparent  soil  thermal  conductiv- 
ity as  a  function  of  moistiu*e  content  and  tempera- 
ture (Campbell  and  others  1994)  and  a  new  water 
content-humidity  relationship  for  soils  (Campbell 
and  others  1993).  It  incorporates  the  advances  in 
theoretical  and  empirical  modeling  of  transport 
properties  over  the  past  20  years  (Campbell  1985). 

Campbell  and  others  (1992,  1995)  model  sensible 
heat  accumulation  as  the  difference  between  the  di- 
vergence of  the  heat  flux  and  the  rate  of  heat  ab- 
sorption per  unit  volume  by  moisture  evaporation. 
This  is  a  fundamental  difference  from  the  approach 
taken  by  Aston  and  Gill.  The  heat  flux  is  modeled 
as  the  product  of  apparent  thermal  conductivity  and 
temperature  gradient,  and  all  soil  components  are 
assumed  to  have  the  same  temperature  at  the  same 
"point"  in  the  continuum.  Liquid  water  is  treated  as 
fi-ozen  in  place,  with  its  concentration  changing  be- 
cause of  evaporation.  This  simplification  is  justified 
on  the  grounds  that  liquid  water  movement  is  too 
slow  to  maintain  equilibrium  with  the  rapid  tem- 
perature changes  that  are  to  be  expected  in  the 
heating  environment  under  a  fire. 

The  flux  density  of  water  vapor  is  calculated  as 
the  negative  of  the  product  of  vapor  conductivity  £md 
the  gradient  of  p,  the  partial  pressure  of  water  va- 
por, all  multiplied  by  the  "Stefan  correction"  P/iP-p), 
where  P  is  total  atmospheric  pressure.  (This  factor 
may  or  may  not  be  appropriate  over  some  range  of 
water  vapor  partial  pressure  but  is  certainly  not 
valid  near  the  boiling  point  of  water,  where  p  =  P 
and  the  equation  for  vapor  flux  becomes  singular. 
See  appendix  A  for  discussion  of  this  factor).  If  we 
assume  that  the  liquid  water  doesn't  move,  and  ne- 
glect the  rate  of  accumulation  of  water  vapor,  with 
respect  to  the  other  terms  in  the  equation,  the  con- 
servation of  mass  applied  to  water  in  the  soil  yields 
an  equation  linking  the  time  rate  of  change  of  liquid 
water  volume  firaction  to  the  divergence  of  the  water 
vapor  flux  density.  The  liquid  water  fraction  is  in 
turn  equated  to  a  compact  empirical  function  of  the 
water  (matric)  potential  that  includes  only  two  pa- 
rameters: the  water  potential  of  ovendry  soil  and  an 
empirical  constant  found  to  be  approximately  six 
times  the  air  dry  moisture  fraction  (Campbell  and 
others  1993).  This  relationship  alone  is  an  impor- 
tant simplifying  advance  in  the  modeling  process. 
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The  apparent  thermal  conductivity  model  is  an 
extension  of  deVries'  (1963)  model  and  is  fully  docu- 
mented for  the  first  time  in  Campbell  and  others 
(1994).  It  is  complicated  but  straightforward,  and  in- 
corporates latent  heat  transport  by  vapor  fl\ix. 

The  vapor  conductivity  model  follows  the  formula- 
tion of  Philip  and  deVries  (1957),  including  the  con- 
catenation of  empirical  factors  multiplying  the  binary 
diffusion  coefficient  of  water  vapor  in  air. 

An  important  contribution  in  this  model  is  its  nu- 
merical integration  scheme,  which  uses  temperature 
and  water  potential  as  independent  variables  in  ex- 
pressions for  the  correction  of  the  errors  in  the  heat 
and  mass  conservation  equations  at  each  node.  The 
integration  scheme  is  fast  and  numerically  stable 
over  a  wide  range  of  test  conditions. 

This  model  represents  an  important  advance  and 
will  be  tested  in  field  appUcations  when  it  is  available. 

Other  models— Three  other  models,  taken  from 
graduate  theses,  are  firmly  in  the  family  of  models 
from  soil  science,  and  the  influence  of  the  pioneering 
works  is  evident. 

Schroeder  (1974)  developed  an  optimized  computer 
simulation  model  for  heat  and  moisture  transfer  in 
soils  for  a  Ph.D.  degree  in  computer  science  at  Texas 
A&M  University.  In  selecting  a  formulation  to  be 
implemented  as  the  computer  simulation,  Schroeder 
summarized  the  equations  describing  moisture  move- 
ment found  in  five  publications.  The  only  model  that 
included  heat  transport  was  the  model  of  Philip  and 
deVries  (1957)  and  deVries  (1958).  The  deVries  (1958) 
model  is  the  one  he  chose  to  implement.  Emphasis 
was  on  the  efficiency  of  the  computer  code  imple- 
menting the  model,  and  not  the  model  itself  There 
was  no  extension,  development,  or  experimental 
verification  of  the  model. 

In  contrast,  the  dissertation  of  Jury  (1973),  docu- 
ments the  writing  and  testing  of  a  computer  simula- 
tion, also  discusses  the  physics  underlying  the  pro- 
cesses modeled,  and  compares  model  predictions 
with  experimental  results.  Cited  as  motivation  for 
investigating  the  workings  of  the  deVries'  model 
(which  is  essentially  the  model  he  used)  was  in  dis- 
agreement between  experiment  and  theory  (Jury 
1973,  p.  4): 

. .  .Several  experiments  have  cast  serious  doubt  on 
the  present  theoretical  models  of  the  coefficients 
used  for  describing  the  thermal  transport  of  water 
in  the  liquid  phase.  None  of  the  models  has  yet  been 
able  to  predict  detailed  behavior  in  the  vapor  phase. 

Jury  presented  a  flowchart  and  FORTRAN  listing 
of  the  program  developed,  but  this  particular  model 
was  not  implemented  for  testing  because  it  was  not 
directed  toward  the  temperature  and  heating  regime 
associated  with  a  fire  environment.  Nonetheless, 


Jury's  thesis  remains  valuable  for  the  physical  con- 
cepts behind  the  mathematical  expressions  in  these 
complicated  models. 

Peter  (1992)  directly  attacked  one  of  the  situations 
that  motivate  the  entire  realm  of  study  here:  heat 
transfer  in  soil  beneath  a  spreading  fire.  Peter  docu- 
ments four  major  areas  of  study: 

1.  Heat  transfer  within  fires  in  a  fuel  bed 

2.  Heat  transfer  within  dry  mineral  soil 

3.  Heat  and  mass  transfer  within  moist  soil 

4.  Heat  and  mass  transfer  within  an  organic  layer 
undergoing  smoldering  combustion  above  mineral 
soil 

Peter's  work  also  includes  coding  of  a  computer  sim- 
ulation. The  basis  for  the  model  is  the  formulation 
of  deVries  (1958),  despite  the  reservations  noted  and 
the  observations  of  deficiencies  made  by  Jury  (1973). 
The  work  leaves  little  room  for  conjecture.  A  work- 
ing personal  computer  version  of  the  simiilation  code 
could  not  be  obtained  for  testing  for  this  study.  The 
numerical  integration  algorithm  it  contains  was  a 
supercomputer-based  utility  routine  that  is  being 
replaced. 

Luikov's  Formulation— No  survey  of  this  field 
would  be  complete  without  acknowledging  the  broad 
contributions  of  the  Russian  scientist,  A.  V.  Luikov, 
head  of  the  Heat  and  Mass  Transfer  Institute  of  the 
Belorus  Academy  of  Sciences.  His  major  works  are 
summarized  in  a  survey  paper  (Luikov  1975),  pub- 
lished posthumously. 

LuikoVs  survey  covers  research  contributions  from 
the  former  Soviet  Union  to  the  field  of  heat  and 
mass  transport  in  porous  media,  by  providing  a  math- 
ematical synopsis  of  models  developed  at  the  Heat 
and  Mass  Transfer  Institute.  While  it  should  have 
worldwide  research  community  backing,  little  of  the 
background  work  is  accessible  in  English  and  the 
svu-vey  is  not  sufficiently  detailed  to  permit  replica- 
tion of  the  computations  exhibited.  Our  unsatisfying 
and  inescapable  conclusion  is  that  the  supporting 
work  must  be  duplicated,  but  that  Luikov's  sum- 
mary guidance  should  be  valuable  in  that  process. 

Models  from  Engineering  and  Geophysics 

Engineering  and  geophysics  literature  (see  cita- 
tions in  Introduction)  shows  contrast  between  ap- 
proaches to  the  formulation  of  models  for  heat  and 
mass  transport  in  porous  media.  Engineering  ap- 
proaches to  modeling  include  Kansa  and  others 
(1977)  who  employed  a  model  for  heat  and  mass 
transport  within  a  moist  wood  element  undergoing 
pyrolysis  to  infer  the  theoretical  influence  of  fuel 
moisture  on  the  rate  of  pyrolysis.  Moallemi  and 
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others  (1993)  modeled  the  effect  of  water  content  on 
the  smoldering  combustion  of  woody  materials.  Kallel 
and  others  (1993)  described  drying  bricks  and  curing 
concrete,  as  had  Eckert  and  Faghri  (1980).  Cheng 
(1978)  and  Garg  and  Kassoy  (1981)  looked  at  geo- 
thermal  energy  in  their  model  development.  Rust  and 
Roberts  (1990)  examined  the  role  of  buoyancy-in- 
duced convective  heat  transfer  in  saturated  soil 
heated  from  above,  through  a  building's  concrete 
slab  floor. 

All  these  authors  start  with  continuum  transport 
equations  for  the  conservation  of  mass  for  each 
phase,  for  the  conservation  of  momentum,  and  for 
the  conservation  of  energy.  The  unique  character 
of  the  porous  medium  is  incorporated  through  the 
constitutive  relationships  that  define  the  continuum. 
Approximations  are  expressed  in  terms  of  influences 
on  the  constitutive  relationships. 

Kaviany  (1991)  and  Cheng  (1978)  are  especially 
instructive  in  terms  of  the  differences  between  ap- 
proaches. Because  the  conservation  laws  expressed 
as  continuum  transport  equations  deal  with  the  con- 
centrations—quantities per  ^xmt  volume— of  con- 
served extensive  properties  of  material  occupying 
the  space  in  question,  it  is  important  to  be  precise 
in  identifying  the  differential  volume  intended  when 
"per  xmit  volume"  is  used.  For  example,  the  term 
"mass  of  water  vapor  per  imit  volume"  could  refer  to 
the  local  density  of  the  water  vapor  in  the  gaseous 
phase,  to  the  local  average  mass  of  water  vapor  per 
unit  interstitial  volume  (some  of  which  could  con- 
tain liquid  water),  or  to  the  local  average  mass  of 
water  vapor  per  unit  volume  of  the  soil  medium  (in- 
cluding interstices  and  whatever  resides  in  them). 
These  obviously  are  not  equivalent  quantities. 

Cheng  introduced  three  mathematical  theorems 
that  relate  rates  of  change  in  space  and  time  of  the 
averages  of  properties  taken  over  a  small— nominally 
differential— volume  of  space,  to  averages  over  the 
same  volume  of  the  rates  of  change  of  the  properties. 
Two  theorems,  one  attributed  to  S.  Whitaker  and 
one  to  W.  G.  Gray,  appeared  in  the  literature  of 
chemical  engineering  research.  These  fundamental 
relationships  are  applied  to  well-known  transport 
equations  applicable  on  the  microscopic  scale  to 
each  phase  (solid,  liquid,  and  gas)  present  in  the  po- 
rous medium.  These  relationships  are  manipulated 
to  derive  macroscopic  volume-average  equations  ap- 
plicable to  the  medium  as  a  continuum,  in  which  ap- 
pear the  properties  of  the  phases  used  to  describe 
them  on  the  microscopic  scale.  Cheng  noted  that  the 
conservation  equations  derived  in  this  manner  can 
also  be  obtained  by  careful  consideration  of  a  differ- 
ential volume  of  the  porous  medium.  Use  of  the  av- 
eraging theorems  keeps  one  aware  of  the  ease  with 
which  flawed  reasoning  can  arise. 


Some  papers  focus  on  the  numerical  aspects  of 
modeling  the  processes  involved  in  heat  and  mass 
trEinsport  in  porous  media.  Some  investigations  are 
substantive  and  even  profound  (Vafai  and  Tien  1989). 
Sometimes  they  are  only  mathematical  exercises 
even  though  they  bear  titles  that  imply  the  physical 
processes  modeled  Eire  themselves  the  subjects  of  in- 
vestigation. For  example,  Thomas  and  others  (1980) 
offler  "a  fully  nonlinear  analysis  of  heat  and  mass 
transfer  problems  in  porous  bodies"  that  is  a  highly 
simplified  application  of  a  formulation  attributed  to 
Luikov  (1966)  to  analyze  the  process  of  kiln-drying 
lumber.  The  conclusion  from  this  paper  is  that  a 
linear  theory  with  fixed  constitutive  parameters  is 
adequate  for  this  application.  Similarly,  Yeh  and 
Luxmoore  (1983)  address  "modeling  moisture  and 
thermal  transport  in  unsaturated  porous  media"  in 
an  investigation  related  to  the  transport  and  dis- 
persal of  waterbome  hazardous  waste.  The  paper 
compares  solutions  of  a  contrived  example  by  two 
different  numerical  methods,  one  of  which  is  called 
"exact,"  the  other  identified  only  by  an  acronym. 

The  effects  of  experimental  apparatus  on  the  out- 
come of  experiments  also  preoccupies  some  authors. 
If  the  concerns  are  valid,  some  of  the  cautions  raised 
should  be  of  interest  to  soil  scientists  even  if  the  range 
of  parameters  is  different.  For  instance,  in  labora- 
tory experiments  involving  the  heating  of  laterally 
confined  moist  soil  samples  from  above  or  below,  the 
region  near  the  lateral  boundary  may  offer  much 
less  resistance  to  the  motion  of  fluid  than  does  the 
bxilk  medium  (Georgiadis  and  Catton  1988).  This 
effect  can  arise  if  the  lateral  boxmdary  is  rigid,  smooth, 
and  impenetrable  by  the  soil  particles.  This  con- 
straint can  result  in  a  looser  packing  of  soil  particles 
over  a  region  of  several  soil  particle  diameters  in  ex- 
tent near  the  boundary,  forming  a  natvu"al  channel 
for  "replacement  air"  to  enter  the  medium  if  buoy- 
ancy causes  the  heated  gas  phase  to  exit  vertically 
from  the  interior  region. 

Differences  in  Approach 

In  the  engineering  and  geophysics  literature,  the 
movement  of  fluids  through  porous  media  is  almost 
always  modeled  by  the  momentum  equation  in  which 
diffusive  mass  transport  is  neglected.  This  equation 
expresses  the  acceleration  of  fluid  in  response  to  the 
forces  acting  upon  it,  where  acceleration  is  set  to 
zero  and  fluid  velocity  is  such  that  the  forces  acting 
on  the  fluid  are  in  balance.  One  such  balance  for  liq- 
uid water  in  saturated  soil  equates  the  water  pres- 
sure gradient  to  the  drag  resistance  force  per  unit 
volume,  which  leads  to  Carey's  law  in  the  limiting 
form  as  velocity  approaches  zero.  In  this  limit,  the 
drag  force  due  to  fluid  flow  over  a  small  particle  is 
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proportional  to  the  velocity  of  the  fluid  rather  than 
to  the  square  of  the  velocity.  In  civil  engineering  lit- 
erature, Darc^s  law  is  extended  to  higher  fluid  ve- 
locities by  multiplying  the  hydrauhc  conductivity  by 
a  factor  involving  the  square  of  the  fluid  velocity,  and 
known  as  the  Forchheimer  correction.  Many  other 
correction  factors,  modifications,  and  extensions  of 
Darc/s  law  all  stem  from  analysis  of  the  momentum 
balance  for  fluids  flowing  through  porous  media.  In 
the  soil  science  literatvire,  Darc/s  law  is  properly  iden- 
tified as  empirical,  but  no  connection  to  fluid  mo- 
mentum conservation  is  made. 

In  the  soil  science  formulation  of  Darc/s  law,  pres- 
sure is  replaced  by  the  "matric  potential"  that  is  es- 
sentially the  Gibbs  Free  Energy  (PhiUp  1975).  This 
may  be  a  convenient  way  to  introduce  the  force  asso- 
ciated with  attraction  between  soil  particles  and  liq- 
uid water  such  as  are  involved  in  capillary  attraction 
and  surface  film  formation.  But  by  doing  so,  tem- 
perature and  liquid  water  content  sensitivities  of  a 
highly  nonlinear  character  and  of  uncertain  depen- 
dence on  the  nature  of  the  soil  medium  are  intro- 
duced where  only  the  well-known  temperature  de- 
pendence of  the  viscosity  of  water  existed  before.  If 
the  force  associated  with  attraction  between  liquid 
water  £ind  soil  matrix  had  been  represented  by  a  sepa- 
rate term  in  the  momentum  equation,  it  would  be 
conceptually  clear,  though  the  final  result  could  math- 
ematically be  the  same.  A  threshold  effect,  in  which 
the  attraction  between  water  and  quartz  solid  van- 
ishes abruptly  near  65  °C  (Churaev  1975),  would  be 
easier  to  quantify  if  separated  from  the  masking  ef- 
fect of  diminishing  water  viscosity  with  increasing 
temperature. 

Another  term  appearing  in  the  momentum  equa- 
tion is  the  body  force  due  to  buoyancy,  which  arises 
fi"om  the  difference  between  the  local  static  pressure 
gradient  and  the  local  fluid  density.  If  the  porous 
medium  is  open  laterally  to  the  ambient  pressure 
field,  then  its  gradient  establishes  the  static  back- 
ground within  which  a  density  defect  results  in  a 
net  force  per  unit  volume  equal  to  the  microscale  av- 
erage density  defect  times  the  acceleration  of  grav- 
ity. Even  if  the  local  soil  environment  is  physically 
isolated  from  the  ambient  pressure  field,  a  buoyancy 
force  can  yet  arise  fi-om  a  nonuniform  temperature 
field,  if  the  density  of  the  fluid  changes  with  tem- 
perature. The  effect  of  buoyancy  is  important  in  geo- 
thermal  processes  and  is  often  manifested  in  Benard 
convection  cells  that  transport  large  volumes  of  liq- 
uid water  (Cheng  1978;  Garg  and  Kassoy  1981)  and 
therefore  large  quantities  of  sensible  heat.  But  by 
inducing  liquid  phase  convective  flow,  heat  transfer 
can  be  significantly  altered  (Rust  and  Roberts  1990). 

A  buoyancy  force  should  also  be  in  evidence  in  the 
motion  of  water  vapor  in  the  heated  soil  under  a 


fire.  Unless  a  one-dimensional  idealization  of  the 
process  is  carried  to  the  extreme  that  air  is  excluded 
fi-om  the  hot  soil  under  the  fire,  so  the  atmosphere's 
ambient  hydrostatic  pressure  field  does  not  penetrate 
into  the  soil,  a  local  density  defect  should  give  rise 
to  a  substantial  vertical  body  force.  This  would  cause 
the  air  and  water  vapor  that  make  up  the  gas  phase 
to  rise.  The  rising  motion  woiild  be  resisted  by  iner- 
tia (including  entrainment  of  stationary  or  slower 
moving  fluid)  and  friction  with  the  soil  medium.  The 
net  movement  of  the  gas  phase  components  would 
include  molecular  diffusion  of  the  constituents  rela- 
tive to  the  moving  bulk  medium. 

If  there  is  a  net  upward  motion  of  the  hot  gas  phase, 
there  may  be  a  replenishing  inflow  of  ambient  air 
that  cannot  be  accommodated  within  the  ideaUza- 
tion  of  a  one-dimensional  model.  The  "replacement 
air"  inflow  might  be  idealized  as  confined  laterally 
to  seams,  fissiures,  or  other  soil  nonuniformities,  with 
lateral  dispersion  at  all  depths  as  demanded  by  the 
continuity  equation  for  the  gas  phase.  Vertical  vari- 
ation of  the  lateral  velocity  field  for  this  flow  could 
be  described  by  the  momentum  equation  integrated 
over  a  lateral  region  bounded  by  edges  along  which 
air  inflow  occurs.  Documenting  this  speculative  dis- 
cussion indicates  that,  while  the  one-dimensional 
idealization  poses  conceptual  problems  to  the  ana- 
lyst, who  would  incorporate  buoyancy  into  the  gas 
phase  momentum  equation,  the  problems  are  not 
seen  as  insurmountable.  By  determining  the  lateral 
scale  of  convective  cells  (Koschmeider  1993)  and  re- 
ducing dimensionality  using  Galerkin  averaging  or 
some  similar  method,  a  self-consistent  and  simpli- 
fied one-dimensional  model  should  be  achieved.  This 
speculation  applies  only  to  models  based  on  the  en- 
gineering pEiradigm,  as  buoyancy  terms  are  neglected 
in  models  in  soil  science.  Vafai  and  Tien  (1989)  re- 
ported their  numerical  investigations  showed  that 
only  at  "high"  pressure  gradients  can  a  one-dimen- 
sional approximation  be  used  for  combined  heat  and 
mass  transport  in  porous  media. 

Often  soil  science  models,  like  engineering  models, 
include  the  effect  of  gravity  on  the  liquid  water  ph£ise. 
The  phase  composition  of  the  interstitial  fluid,  in 
the  region  where  liquid  water  exists,  can  vary  from 
all  gas  to  all  liquid,  making  the  modeling  of  this  re- 
gion one  of  the  most  intricate  and  challenging.  It 
should  be  possible  to  represent  the  interstitial  fluid 
as  a  single  continuum  consisting  of  two  intermingled 
phases  whose  relative  abundance  varies  spatially 
and  temporally  (Vafai  and  Tien  1989).  Instead,  sep- 
arate models  for  the  gas  and  liquid  phases  are  de- 
veloped and  combined  with  the  dismissal  of  rela- 
tively small  terms.  A  comprehensive  model  of  the 
two-phase  regime  may  reside  in  the  petroleum  engi- 
neering literatm-e. 
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The  movement  of  water  vapor  in  soil  is  modeled  as 
molecular  diffusion  in  all  the  formulations  from  the 
soil  science  literatvu-e.  Molecular  diffusion  is  a  process 
of  mixing  in  which  the  velocity  of  one  species  of  mol- 
ecule relative  to  another  through  which  it  is  moving 
is  proportional  to  the  gradient  of  the  molecular  con- 
centration of  the  diffusing  species.  These  models  are 
almost  always  formulated  in  terms  of  the  gradient 
of  water  vapor  pressure  or  relative  humidity  and 
usually  include  a  nonlinear  enhancement  factor— 
the  Stefan  correction— that  may  not  be  appropriate 
(see  appendix  A).  The  vapor  pressure  is  in  turn  re- 
lated to  the  local  temperature  so  the  concentration 
gradient  is  determined  by  the  local  temperature 
gradient,  and  the  process  of  heat  transfer  by  vapor 
transport  is  seen  phenomenologically  to  mimic  ther- 
mal conduction.  The  complexity  of  this  connection 
prevents  its  expUcit  modeling,  but  the  analysis  serves 
as  a  guide  to  definition  and  interpretation  of  the  ap- 
propriate empiricism. 

Expressing  the  conservation  of  energy  by  equating 
the  accumulation  of  heat  per  unit  volume  to  the  neg- 
ative of  the  divergence  of  the  rate  of  transport  of  heat 
via  mass  motion  and  bulk  medium  thermal  conductiv- 
ity closes  the  equation  set.  It  also  strongly  couples 
the  models  for  heat  and  mass  transport.  As  the  wa- 
ter vapor  diffusion  model  is  a  strong  function  of  tem- 
perature, movement  of  water  vapor  in  response  to  tem- 
perature gradients  theoretically  results  in  strong 
temperature  dependence  of  the  apparent  thermal  con- 
ductivity of  the  air  that  fills  the  interstices  between 
soil  particles  not  filled  with  liquid  water.  Because 
the  latent  heat  of  vaporization  of  water  multiplied 
by  the  mass  flux  of  water  vapor  represents  the  trans- 
port of  energy,  while,  as  noted,  diffusion  of  water  va- 
por is  modeled  in  part  as  nonlinearly  proportional  to 
temperature  gradient. 

Because  of  the  uncertainty  and  complexity  of  the 
dependence  of  heat  transport  upon  water  vapor  move- 
ment, any  model  must  rely  to  some  degree  upon  an 
empirically  determined  apparent  thermal  conductiv- 
ity that  depends  on  temperature,  water  content,  and 
soil  characteristics.  This  procedure  might  predict 
transient  temperature  fields  with  high  fidelity,  yet 
fail  to  describe  transient  moisture  distributions  ac- 
curately. The  heat  transfer  boundary  conditions  de- 
termine the  temperature  field  from  which  the  mois- 
ture transport  can  be  inferred.  Though  oversimplified, 
the  characterization  that  the  soil  science  models  de- 
rive the  moisture  distribution  from  the  temperature 
field  is  not  strongly  at  variance  with  the  literature 
reviewed. 

Engineering  and  soil  science  differences  are  evi- 
dent in  the  treatment  of  heat  transfer.  The  formula- 
tion based  on  continuum  mechanics  distinguishes 


between  energy  transported  by  mass  motion  and 
heat  transferred  through  the  constituents  of  the  me- 
dium, Temperatiire  fields  for  each  of  the  components 
(solid  particles,  liquid  water,  interstitial  gas  com- 
posed of  air  and  water  vapor)  would  be  determined 
from  the  heat  trgmsferred  to  and  from  each  phase 
and  the  transport  of  heat  by  its  mass  motion.  The 
temperature  of  each  of  the  components  need  not 
necessarily  be  the  same  everywhere.  Temperature 
differences  could  result  in  significant  heat  transfer 
by  vapor  transport  through  and  condensation  on 
cooler  soil  particles  relatively  far  from  the  point  at 
which  the  water  evaporated.  Likewise,  superheating 
of  vapor  rising  through  hot,  dry  overburdened  soil 
could  reduce  heat  flux  to  the  moist  soil  below. 

A  model  constructed  from  a  more  general  perspec- 
tive could  offer  some  advantages.  It  should  collapse, 
in  the  proper  limits,  to  a  diffusion  driven  (soil  sci- 
ence type)  model  or  to  an  engineering  model  (with 
diffusive  terms  neglected)  covering  the  same  phe- 
nomenology. It  should  contain  within  it  the  compet- 
ing terms  that  are  dismissed,  allowing  quantitative 
definition  of  the  ranges  of  validity  of  the  specialized 
limiting  forms. 


Computer  Programs  Exercised 

The  Montana  State  University  contingent  converted 
two  models  to  FORTRAN  suitable  for  application  on 
personal  computers  to  facilitate  testing  of  the  mod- 
els for  their  accuracy  in  predicting  results  of  experi- 
ments at  the  Intermountain  Fire  Sciences  Labora- 
tory. We  present  results  from  testing  two  models, 
Aston  and  Gill  (1976)  and  Campbell  and  others  (1992). 

The  model  of  Aston  and  Gill  (1976)  as  supplied 
by  the  Fire  Lab  to  the  University  study  team  was 
transcribed  to  a  dialect  of  FORTRAN  acceptable  to 
Microsoft  FORTRAN  5.0.  The  model  was  exercised 
to  create  temperature  and  soil  moisture  field  histo- 
ries (examples  of  the  output  in  figs.  1  and  2).  The 
aberrant  behavior  of  the  temperature  histories  at 
the  shallower  depths  is  not  unique  to  the  new  ver- 
sion of  the  model.  The  same  behavior  had  been  noted 
by  Fire  Lab  researchers  in  exercising  a  version  of 
this  model  adapted  to  the  Data  General  at  the  lab. 

Further  testing,  relying  on  soil  constitutive  pa- 
rameter tables  generated  using  Campbell's  (1985) 
predictive  models,  confirmed  that  the  numerical  in- 
tegration scheme  used  in  the  Aston-Gill  model  be- 
comes unstable  when  a  simulated  soil  layer  becomes 
nearly  dry  (Albini  and  Amin  1994).  The  equations 
are  remarkably  "stiff,"  a  characteristic  that  often 
forces  use  of  special  nxmierical  techniques  (Press  and 
others  1986).  This  feature  was  noted  by  the  model 
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Figure  1 — Temperature  histories  at  various  deptins  predicted 
by  Aston-Gill  model.  Depth  1  is  0.5  cm  below  the  surface, 
depth  2  at  1 .5  cm,  and  depth  3  at  2.5  cm.  The  surface  temp- 
erature is  a  rising  exponential  followed  by  a  decaying  expo- 
nential in  time,  much  like  the  temperature  history  at  depth  1. 
Note  the  aberrant  behavior  between  2,000  and  3,000  s. 
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Figure  2 — Temperature  histories  at  greater  depths,  as 
predicted  by  the  Aston-Gill  model  for  the  conditions  of  figure  1 . 
Depth  4  is  3.5  cm  beneath  the  surface,  depth  5  at  4.5  cm,  and 
depth  6  at  5.5  cm.  These  histories  do  not  misbehave  as  badly 
as  do  those  shown  in  figure  1 . 


builders,  but  they  used  a  second  order  Runge-Kutta 
integration  scheme  that  is  not  robust  against  the 
severely  nonlinear  dependencies  arising  in  many 
cases  of  soil  heating  by  fire. 

The  model  of  Campbell  and  others  (1992)  was 
provided  to  the  University  study  team  in  Turbo 
Pascal  source  code.  It  was  translated  into  Microsoft 
FORTRAN  5.0,  compiled,  and  exercised  (example  out- 
put in  figs.  3  to  6).  These  results  simulate  experiments 
in  which  the  upper  surface  of  a  soil  sample  was  ex- 
posed to  a  powerful  radiant  heater.  The  upper  sur- 
face temperature  histories  were  brought  into  agree- 
ment by  adjusting  the  parameters  that  describe  the 
upper  surface  heat  input  rate.  With  these  param- 
eters fixed,  predictions  of  thermal  and  moisture  re- 
sponse were  made  for  a  wide  variety  of  soils  under 
various  initial  moisture  conditions.  Measurements 
of  the  temperature  at  various  depths  showed  good 
agreement  with  the  model's  predictions  (Campbell 
and  others  1992,  1995). 

Campbell  and  others  (1992,  1995)  compare  numer- 
ous experimental  and  theoretical  results;  in  almost 
every  case,  the  model  predicts  at  least  the  high  tem- 
perature field  history  with  good  accuracy,  but  yields 
relatively  large  errors  in  prediction  of  moisture  con- 
tent histories.  A  later  variant  of  this  model  (Campbell 


and  others  1995)  should  improve  the  accuracy  of  mois- 
ture content  predictions.  The  measurement  of  local 
moisture  concentration  in  soils,  even  in  the  labora- 
tory, is  difficult  without  disturbing  the  experimental 
conditions.  The  temporal  and  spatial  resolutions  of 
nonintrusive  instruments  are  inadequate,  and  high 
resolution  methods  are  intrusive. 

The  Aston  and  Gill  (1976)  model  is  not  in  a  usable 
state  at  this  time.  The  Campbell  and  others  (1992) 
model  is  now  available  in  two  source  languages,  is 
operable  on  personal  computers,  and  is  fast  and  easy 
to  use  in  either  variant.  It  seems  to  perform  well  in 
describing  the  transient  temperature  fields  found  in 
laboratory  experiments.  Although  its  prediction  of 
moisture  movement  is  not  well  validated  by  labora- 
tory measurements,  it  may  satisfy  all  requirements 
for  a  general  heat  transfer  model  for  fire  effects  pre- 
diction. It  clearly  should  be  tested  for  this  purpose. 


Opportunities  for  Model  Improvement 

Opportunities  for  possible  improvement  in  the 
modeling  of  heat  and  mass  transport  in  soils  heated 
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Figure  3 — Temperature  histories  predicted  by  the  Campbell 
model  for  soil  surface  and  two  depths  beneath  the  surface. 
Radiant  heating  by  a  high  powered  lamp  in  a  laboratory 
experiment  is  simulated.  The  soil  is  air  dry  Boulder  Creek  silt 
loam  at  an  initial  moisture  content  of  6  percent. 
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Figure  4 — Moisture  content  histories  corresponding  to  the 
temperature  histories  displayed  in  figure  3,  as  predicted  by 
Campbell's  model.  Temporary  rise  in  moisture  content  at 
25  mm  depth  indicates  an  initial  vapor  migration  downward. 
Zero  moisture  content  refers  to  105  °C  ovendry  condition,  so 
negative  moisture  predictions  are  physically  meaningful. 
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Figure  5 — ^Temperature  histories  predicted  by  Campbell's 
model  for  the  same  soil  and  heating  conditions  of  figure  3, 
but  with  initial  soil  moisture  of  22  percent.  Note  the  higher 
temperatures  at  65  mm  depth,  presumably  due  to  higher 
thermal  conductivity  in  the  wetter  soil. 


by  fires  must  be  viewed  as  hypothetical.  They  are  of- 
fered here  as  suggestions  for  the  creators  of  the  next 
generation  of  models.  They  fall  into  two  categories: 

1.  Identification  of  phenomena  left  out  of  models 
that  might  prove  to  be  important. 

2.  Identification  of  complicating  features  included 
that  might  prove  to  be  unnecessary. 


Possible  Important  Omissions 

Models  developed  for  engineering  applications  al- 
most always  ignore  diffusive  transport  of  water  as 
vapor  or  as  liquid,  and  so  lack  representation  of  this 
phenomenology.  Such  models  can  only  be  accurate 
in  the  limiting  situation  that  diffusive  mass  motions 
can  be  neglected.  Since  the  applicability  of  the  limit 
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Figure  6 — Moisture  histories  corresponding  to  the  temperature 
histories  shown  in  figure  5.  Transient  increase  in  moisture 
content  is  predicted  for  25  mm  depth.  Zero  moisture  content 
refers  to  105  °C  ovendry  condition,  so  negative  moisture 
predictions  are  physically  meaningful. 


cannot  be  inferred  from  the  content  of  the  models, 
they  are  fundamentally  deficient  for  this  application. 
Therefore,  we  address  omissions  in  the  formulations 
of  models  with  soil  science  ancestry. 

Omission  of  the  momentum  equation  from  the  prim- 
itive set  appears  to  be  unjustified.  Clearly,  if  heat  is 
supplied  rapidly  to  very  moist  soil,  the  soil  mediimi 
eventually  will  emit  steam  at  a  rate  that  cannot  be 
carried  away  by  the  process  of  vapor  diffusion.  A 
pressure  gradient  would  develop,  negative  upward, 
that  would  assist  the  (omitted)  buoyancy  force  in  ex- 
pelling the  steam  in  a  massflow  process.  This  flow 
would  inhibit  the  transfer  of  heat,  and  a  pseudosteady 
state  is  conceivably  achievable. 

The  problem  of  predicting  the  transient  movement 
of  the  phase-change  boundary  (the  classical  "Stefan 
problem")  is  modified  by  the  heat  transported  by  the 
vapor  phase  massflow  counter  to  the  flow  of  heat 
from  fire  to  moist  soil.  It  could  even  require  an  eigen- 
value solution.  This  would  be  a  most  interesting 
theoretical  development  and  could  have  quite  im- 
portant fire  effects  implications  in  the  field. 

The  assumption  that  all  components  of  the  soil 
medium  have  the  same  temperature  at  each  point  has 
strong  heuristic  appeal.  The  assumption  may  have 
arisen  from  the  piecemeal  assembly  of  an  energy 
transport  equation  starting  with  Fourier  heat  con- 
duction and  an  apparent  thermal  conductivity  for 
the  medium.  The  variations  in  this  equation  among 
model  authors  arises  from  the  primitive  equation  set 
chosen,  which  is  not  a  fundamental  one.  A  continuum 
mechanics  approach,  with  multiple  fluid  components 
occupying  the  interstices  of  a  solid  matrix,  would 
readily  allow  separate  temperature  distributions  for 
each  of  the  components.  The  heat  transfer  rates  be- 
tween components  would  have  to  be  modeled.  Such  a 
model  should  predict  the  conditions  under  which  the 
assumption  is  valid  that  all  components  locally  have 
the  same  temperature. 


We  did  not  discover  a  treatment  of  lateral  non- 
homogeneity  from  soil  science  sources,  nor  did  we 
find  an  example  of  a  systematic  reduction  of  dimen- 
sionality by  starting  with  a  three-dimensional  model 
and  applying  averaging  to  reduce  the  model  to  one 
dimension.  This  is  true  despite  watershed  works  of 
Philip  and  deVries  (1957)  and  of  deVries  (1958)  that 
formulated  the  process  descriptions  in  three-dimen- 
sional form.  This  could  contribute  to  the  difficulties 
in  predicting  moisture  movement  that  have  beset 
these  models  (Campbell  and  others  1992;  Jury  1974). 
This  implication  arises  because  a  one-dimensional 
model  imposes  the  decision  as  to  the  direction  of  va- 
por movement  to  be  based  on  possibly  competing 
responses  to  thermal  and  moisture  gradients.  This 
allows  the  possibility  of  a  zero-velocity  front  that 
marks  a  boundary  between  vapor  moving  in  the  di- 
rection of  heat  flow  (aiding  heat  transfer)  and  mov- 
ing in  the  direction  covmter  to  heat  flow.  We  have  not 
discovered  a  thorough  discussion  of  this  phenom- 
enology. A  model  that  countenances  lateral  move- 
ment of  water  vapor  and  interstitial  air  might  alle- 
viate such  complications. 

Finally,  we  have  seen  no  accommodation  or  repu- 
diation of  the  effect  documented  by  Churaev  (1975) 
that  the  wetting  attraction  of  liquid  water  to  the  in- 
terior surfaces  of  quartz  capillary  tubes  rapidly  goes 
to  zero  with  increasing  temperature  near  65  °C.  If 
this  phenomenon  occurs,  does  it  occur  in  sandy  soils 
also?  If  not,  why  not?  If  it  occurs,  there  should  be  a 
different  dependence  of  water  matric  potential  on 
moisture  content  at  high  temperature  than  at  low. 
This  potentially  important  phenomenon  appears  to 
have  been  ignored.  Either  that,  or  we  have  failed  to 
unearth  the  definitive  dismissal  of  the  importance 
of  this  phenomenon  to  moisture  movement  in  heated 
soils.  This  phenomenon  needs  to  be  investigated 
further. 
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Possible  Simplifications 

Models  from  soil  science  ancestry  are  complex  and 
could  be  simplified  without  seriously  degrading  their 
capabilities  to  predict  heat  and  mass  transport  in 
soils  heated  by  fires.  The  modifications  necessary  to 
render  engineering  models  applicable  to  fire-driven 
heat  and  mass  transport  in  soils  amount  to  the  ad- 
dition of  contributions  to  mass  transport. 

The  apparent  thermal  conductivity  used  to  model 
heat  transport  in  soils  employed  by  deVries  (1958) 
is  more  clearly  rationalized  by  Campbell  and  others 
(1994).  In  the  earlier  formulation,  latent  heat  trans- 
port by  vapor  movement  is  modeled  by  invoking  the 
condition  that  vapor  movement  is  controlled  by  tem- 
perature dependence  of  the  rate  of  diffusion  of  water 
vapor.  Using  this  formulation  makes  the  apparent 
thermal  conductivity  a  strongly  nonlinear  function 
of  moisture  content  and  of  temperature  and  requires 
the  use  of  a  source  or  sink  term  in  the  Fourier  equa- 
tion employing  the  apparent  conductivity.  The  pre- 
dicted contribution  of  vapor  movement  to  apparent 
thermal  conductivity  seems  to  require  an  empirical 
adjustment  that  increases  it  by  an  order  of  magnitude. 
Under  the  high  heating  rates  and  temperatures  of 
fire-driven  heat  transfer,  it  is  conceivable  that  a  for- 
mulation based  on  a  more  primitive  equation  set 
would  be  simpler.  The  movement  and  temperature 
of  each  of  the  phases  in  the  bulk  medium  would  be 
modeled  separately,  taking  into  account  heat  trans- 
ferred between  phases.  Limiting  case  approximations 
might  be  more  transpeu-ent  under  such  a  formulation. 

The  descriptions  of  vapor  movement  used  in  the 
models  based  on  soil  science  rest  upon  Pick's  law  for 
binary  molecular  diffusion.  This  law  states  that  the 
net  flux  of  molecules  (number  of  molecules  per  unit 
area  per  unit  time)  is  proportional  to  the  negative 
gradient  of  the  concentration  of  such  molecules.  Ap- 
plied to  the  interstitial  gas  phase  treated  as  a  binary 
mixture  or  air  and  water  vapor,  yields  a  macroscopic 
phenomenological  equation  that  has  two  multiplica- 
tive empirical  factors:  the  "tortuosity"  correction, 
usually  taken  to  be  about  two-thirds,  and  an  "en- 
hancement" factor  that  ranges  between  5  and  10.  In 
addition,  the  "Stefan  correction"  (see  appendix  A)  is 
applied,  which  increases  the  calculated  flux  of  water 
vapor  molecules  by  a  factor  proportional  to  the  ratio 
of  the  ambient  pressure  to  the  partial  pressure  of 
air.  When  the  local  temperature  approaches  the 
boiling  point  of  water,  this  factor  becomes  infinite  if 
there  is  any  liquid  water  present,  so  it  must  be  lim- 
ited to  some  maximum  value,  usually  less  than  10. 

In  the  engineering  literature,  water  vapor  move- 
ment in  porous  media  is  treated  as  part  of  the  gas 
phase.  The  temperature  of  and  the  local  heat  trans- 
fer rate  to  the  liquid  water  component  and  the  local 


concentration  of  water  vapor  in  the  gas  phase  are 
modeled  as  influencing  the  local  rate  of  evaporation. 
The  gas  phase  moves  through  the  interstitial  pores 
of  the  soil  matrix  that  are  not  blocked  by  liquid  wa- 
ter at  a  rate  fixed  by  the  local  pressiire  gradient,  gas 
friction  with  the  stationary  matrix,  and  body  forces 
acting  on  it.  This  model  results  in  a  high-gas  phase 
velocity  on  the  microscopic  scale,  relative  to  the  mass 
motion  velocity  on  the  macroscopic  scale.  Tempera- 
ture sensitivity  of  the  viscosity  of  the  gas  phase  and 
coupling  of  the  heat  transfer  and  mass  motion  equa- 
tions render  this  model  fully  as  nonlinear  and  as 
complicated  as  the  soil  science  model,  if  it  is  reduced 
by  combining  equations  to  derive  an  equivalent  to 
the  "vapor  conductivity"  used  in  soil  science  models. 
Because  each  step  in  the  process  can  be  described 
separately,  the  model  is  conceptually  simpler  and 
its  failure  points  more  readily  identifiable  when 
comparing  experimental  data  to  predictions.  The 
diffusion  of  water  vapor  in  this  model  would  only 
serve  to  disperse  the  water  vapor  through  the  gas 
phase  and  would  not  be  modeled  explicitly  but  would 
be  assumed  to  have  resulted  in  local  lateral  unifor- 
mity. Contributions  to  mass  motion  from  diffusive 
water  vapor  movement  are  neglected  in  these  for- 
mulations, so  the  Stefan  correction  would  play  no 
role.  Neglect  of  vapor  diffusion  would  mean  that  the 
model  should  not  collapse  to  one  derived  from  soil 
science  in  the  limit  of  small  thermal  gradients  and 
heat  transfer  levels. 

Some  soU  science  models  account  for  the  movement 
of  liquid  water  and  the  transport  of  sensible  heat 
thereby  (Aston  and  Gill  1976;  Peter  1992),  while 
Campbell  and  others  (1992,  1995)  model  liquid  wa- 
ter as  being  "frozen  in  place"  in  fire-heated  soils. 
This  simplification  seems  justified  given  that  the 
flux  of  liquid  water  in  soil  is  governed  by  the  prod- 
uct of  hydraulic  conductivity  and  gradient  of  water 
potential.  The  highly  nonlinear  behavior  of  the 
hydraulic  conductivity  with  soil  moisture  content 
introduces  mathematical  sensitivity  near  the  dry 
condition  that  may  require  special  attention  in  mod- 
eling. Yet  the  movement  of  liquid  water,  even  at  ex- 
treme temperature  and  moisture  gradients,  is  slow 
relative  to  the  movement  of  heat  (as  measured  by 
the  rate  of  advance  of  an  isotherm  at  less  than  the 
boiling  point),  and  its  motion  transports  little  heat. 
Retention  of  this  assumption  would  be  helpful  in 
assembling  a  heat  and  mass  transport  model  from  the 
continuum  mechanics  approach.  Overlooking  some 
important  phenomenology  associated  with  a  temper- 
ature threshold  for  the  surface  wetting  attraction 
between  liquid  water  and  soil  minerals  (Churaev 
1975)  would  justify  retaining  a  model  for  liquid  wa- 
ter movement  in  the  circumstances  that  the  model 
is  intended  to  be  used. 
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Conclusions 

Our  literature  survey  and  study  of  models  for  heat 
and  mass  transport  in  soils  exposed  to  heating  rate 
and  temperature  regimes  expected  under  wildland 
fires  has  led  us  to  conclude: 

•  The  soil  science  field  has  contributed  the  only 
useful  models  for  phenomenology  readily  avail- 
able to  investigators  of  wildland  fire  effects. 

•  The  model  of  Campbell  and  others  (1992,  1995) 
seems  to  perform  well  in  predicting  temperature 
histories,  at  various  depths,  of  soils  heated  on 
their  upper  surfaces  at  rates,  and  to  tempera- 
tures expected  under  wildland  fires.  It  does  not 
seem  to  predict  soil  moisture  content  histories 
nearly  as  well  as  it  does  temperature  histories, 
at  least  in  the  early  variant  examined  here. 

•  An  alternative  model  can  be  described  as  a  con- 
tinuum mechanics  approach  usually  used  within 
the  engineering  disciplines.  This  uses  a  more 
fundamental  set  of  primitive  equations  and  could 
more  readily  allow  for  nonequilibrium  interac- 
tions between  the  constituents  of  the  soil  medium. 
Such  a  model  should  be  modified  to  include  dif- 
fusive mass  motion  of  water  vapor  and  liquid 
water  as  well. 

•  We  have  identified  the  omission,  in  soil  science 
models,  of  several  important  features.  Conversely, 
some  complicating  aspects  of  such  models  might 
contribute  only  negligibly  to  heat  or  mass  trans- 
port. Their  dismissal  might  simplify  the  model- 
ing with  little  sacrifice  in  fidelity. 
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Appendix  A:  The  Stefan  Correction  for  Confined  One-Dimensional  Diffusion 


Consider  the  steady  state  diffusion  of  water  vapor  through  air  under  the 
following  constraints  and  assumptions: 

1.  The  flow  takes  place  in  a  small  diameter  horizontal  tube,  so  motion  is 
confined  to  one  direction  and  gravity  plays  no  role. 

2.  The  left  end  of  the  tube,  at  jc  =  0,  is  immersed  in  liquid  water  so  there  can 
be  no  movement  of  air  through  that  boundary. 

3.  The  right  side  of  the  tube  is  open  to  the  ambient  atmosphere,  so  the  fluid 
in  the  tube  supports  no  pressure  gradient  along  its  length  under  steady  state 
conditions. 

4.  Ambient  relative  humidity  is  less  than  100  percent,  so  water  evaporates 
fi-om  the  liquid  phase  surface  and  escapes  into  the  atmosphere. 

5.  Both  air  and  water  vapor  can  be  ideaUzed  as  perfect  gases,  and  the  pro- 
cesses within  the  tube  take  place  at  ambient  temperature. 

If  water  vapor  is  to  diffuse  along  the  length  of  the  tube,  moving  from  left  to 
right,  then  there  is  a  gradient  in  the  concentration  of  water  molecules,  increas- 
ing from  right  to  left.  There  must  also  be  a  countervailing  gradient  in  the  con- 
centration of  air  molecules,  so  air  diffuses  from  right  to  left.  Since  air  cannot 
escape  the  tube  at  its  left  end,  unlike  the  water  vapor  that  escapes  at  its  right 
end,  there  must  be  a  general  flow  of  the  air-water  vapor  mixture  from  left  to 
right,  keeping  the  net  movement  of  air  at  zero.  This  general  flow  acts  to  in- 
crease the  transport  speed  of  the  water  vapor,  increasing  the  "effective  diffu- 
sion coefficient"  of  the  water  vapor  through  air.  Stefan  (1874)  first  documented 
this  phenomenon.  The  following  derivation  follows  Jakob  (1956). 

Nomenclature: 
C    concentration,  mol/volume 
D   bimolecular  diffusion  coefficient,  lengths/time 
M   molecular  weight,  mass/mol 
m"  mass  flux,  mass/area  x  time 
P    (partial)  pressure,  energy/volume 
R    perfect  gas  constant,  energy/mol  x  degree 
T    (absolute)  ambient  temperature 
U   convection  velocity,  positive  left  to  right,  length/time 
X    distance  along  tube,  positive  left  to  right 

Subscripts 

0  total  (ambient) 

1  refers  to  air 

2  refers  to  water  vapor 

1.  +  P2  =  Po  Constant  pressure 

2.  Ci  =  Pi//?  To  Air  is  perfect  gas 

3.  C2  =  P'2!R  Tq  Water  vapor  is  perfect  gas 

4.  mi"  =  MiiCiU  -  D  dCildx)  =  0  No  net  air  flux 

5.  7712"  =  M2{C2U  -  D  dC2ldx)  Net  water  vapor  flux 

6.  PiU  =  DdPi/dx  4.  and  2. 

7.  P2  f/  =  D  dP2/dx  +  RTo  1712  IM2        5.  and  3. 

8.  PqU  =  RTq  m2lM2  6.  +  7.  and  1. 
Solve  8.  for  U,  to  be  used  in  5.  Solve  resultant  for  1x12  to  find: 

ma"  =  -  [Po/(Po  -  P2)]  M2  D  dCz/dx. 

The  factor  in  square  brackets  is  the  "Stefan  correction"  that  increases  the 
effective  value  of  the  diffusion  coefficient  D.  The  equation  breaks  down  as  the 
value  of  P2  approaches  Pq.  The  flaw  in  reasoning  is  that  the  assumed  steady 
state  will  not  be  achieved. 


16 


Albini,  Frank;  Amin,  M.  Ruhul;  Hungerford,  Roger  D.;  Frandsen,  William  H.;  Ryan, 
Kevin  C.  1996.  Models  for  fire-driven  heat  and  moisture  transport  in  soils.  Gen.  Tech. 
Rep.  INT-GTR-335.  Ogden,  UT:  U.S.  Department  of  Agriculture,  Forest  Service, 
Intermountain  Research  Station.  16  p. 

A  survey  was  conducted  of  predictive  models  for  heat  and  mass  transport  within  soils 
exposed  to  the  heating  rates  and  temperature  regimes  under  wildland  fires.  Two  models 
trace  their  ancestry  to  soil  science,  and  other  models  for  heat  and  mass  transport  in 
porous  media  come  from  engineering  disciplines  and  geophysics.  The  approaches 
underlying  the  development  of  the  models  were  contrasted.  Opportunities  for  improving 
the  rotsustness  and  fidelity  of  currently  available  models  were  sought  and  a  few  potentially 
important  omissions  were  identified,  along  with  some  complicating  features  that  might  be 
ignored  without  significant  loss  of  realism. 
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